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I make the first estimates of maximum elastic quadrupole deformations sustainable by alternatives 
to conventional neutron stars. Solid strange quark stars might sustain maximum ellipticities (dimen¬ 
sionless quadrupoles) up to a few times 10“^ rather than a few times 10“^ for conventional neutron 
stars, and hybrid quark-baryon or meson-condensate stars might sustain up to 10~®. Most of the 
difference is due to the shear modulus, which can be up to 10^^ erg/cm® rather than 10^° erg/cm^ 
in the inner crust of a conventional neutron star. Maximum solid strange star ellipticities are com¬ 
parable to upper limits obtained for several known pulsars in a recent gravitational-wave search by 
LIGO. Maximum ellipticities of the more robust hybrid model will be detectable by LIGO at initial 
design sensitivity. A large shear modulus also strengthens the case for starquakes as an explanation 
for frequent pulsar glitches. 
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The LIGO Science Collaboration recently used data 
from the second LIGO science run (S2) to set upper lim¬ 
its on gravitational wave emission from 28 known pulsars, 
9 of which have no competing upper limit from radio ob¬ 
servations Pi . For those 9 pulsars the best S2 upper lim¬ 
its on neutron star ellipticity [a dimensionless quadrupole 
moment, see Eq. ©] are a few times 10 With better 
data now being analyzed, LIGO will soon be sensitive 
to ellipticities of 10“® or less. This raises the question 
of when a detection might be possible, or when enough 
non-detections (upper limits) begin to confront some the¬ 
oretical models of dense matter. The answer depends on 
the maximum ellipticities sustainable in those models. I 
make the first estimates of maximum ellipticities for sev¬ 
eral exotic matter models and find that the LIGO S2 
search was already sensitive to the upper end of the the¬ 
oretical range. LIGO observational results are becoming 
astrophysically interesting years sooner than previously 
expected. 

The maximum elastic deformation sustainable by a 
neutron star has been addressed several times in the past 
few decades—see P, and references therein. A conven¬ 
tional neutron star consists of a liquid nuclear-matter 
core covered by a thin solid crust, which is responsible for 
the deformation and whose microphysics can be extrap¬ 
olated conservatively from laboratory nuclear physics. 
More exotic models of compact stars have been proposed, 
some including large solid cores (see P for a summary), 
but the maximum deformation has not been quantita¬ 
tively addressed. Historically the problem was of inter¬ 
est first in relation to the “glitch” phenomenon in pul¬ 
sars, which was believed to be related to starquakes 0. 
However, the total elastic energy stored in a maximally 
strained crust is far too low to explain the strength and 
frequency of the glitches of the pulsar Vela X-1 P. Oc¬ 
casionally works on exotic compact stars have mentioned 
that solid cores might revive the starquake glitch mech¬ 
anism, but without estimating numbers. 


In this Letter, I estimate maximum elastic deforma¬ 
tions sustainable by exotic alternatives to neutron stars 
and work out the implications for gravitational-wave 
emission and pulsar glitches. Of the models extant 
in the literature, solid strange stars allow the largest 
ellipticities—up to 10^ times those of neutron stars— 
although this model is highly speculative. Hybrid quark- 
baryon stars and stars with charged meson condensates, 
both based on more robust theories, might allow elliptic¬ 
ities up to a few times 10^ more than those of conven¬ 
tional neutron stars. This makes detectable gravitational 
wave emission a prospect for initial LIGO rather than ad¬ 
vanced LIGO and makes the starquake model of glitches 
viable again. 

There are several sources of uncertainty in such esti¬ 
mates. The largest is the matter model itself—maximum 
ellipticities vary by 10^ between conventional neutron 
stars and solid strange stars. The second largest is the 
breaking strain. I quote fiducial numbers for a breaking 
strain of 10“^, which is near the maximum for terrestrial 
alloys and may be favored by observations of neutron 
stars in low-mass x-ray binaries, but the breaking strain 
could be lower by 10^-10^ 0. For hybrid and condensate 
stars, charge screening might bring the maximum ellip¬ 
ticity down to that of a neutron star. These uncertainties 
justify making several approximations which simplify the 
calculations at a cost of introducing relatively small er¬ 
rors as in P. Relativistic gravity and rotational effects 
can change the density profile of a star by tens of percent; 
but they cancel to some extent and are smaller than the 
effect of varying the star’s mass a few percent P, and 
thus I neglect them. Because of the high Fermi energies 
involved, finite temperature plays a negligible role in de¬ 
termining the maximum ellipticity. In the strange and 
hybrid stars, a normal solid crust is believed to be still 
present, but its contribution to the ellipticity is a few 
percent correction to that of the core. I quote maximum 
ellipticities including the maximum 200% contribution 
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from the self-gravity of the deformation Hi , but that 
could go down by a factor of 2. Further calculations, 
details, and uncertainties will be presented elsewhere Q- 
Neutron stars. —Reference Q computes in its Eq. (69) 
a maximum m = 2 quadrupole moment for a neutron 
star using a chemically detailed model of the crust. Cor¬ 
recting the definition of shear modulus [^, it reads 



where CTmax is the breaking strain of the crust. 

The quadrupole o can be converted to the elliptici^ 
e = {Ixx — lyy)/Izz used in gravitational wave papers [l|: 

e=^J%^T|lhQ22|Izz, ( 2 ) 


where the z-axis is the rotation axis and lab is the mo¬ 
ment of inertia tensor. For conventional neutron stars, 
Bejger and Haensel find that the approximation 
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is accurate to a few percent for a variety of equations of 
state. Thus we can write the maximum ellipticity of a 
conventional neutron star as 



For the fiducial values of mass, radius, and breaking 
strain, Cmax is 2 x 10“^ (6 x 10“^ with self-gravity). 

The generalization of Eq. o to arbitrary equations of 
state can be obtained by combining Eqs. (67) and (64) 
of Ref. 0 as 
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where fi is the shear modulus (nonzero only in the solid 
part of the star), g is the local gravitational acceleration, 
and U = 2 + dliig/dhir. The two bounding cases are in¬ 
compressible matter and infinitely compressible matter 
(a point mass). Note that the latter is equivalent to a 
conventional neutron star, where the mass of the crust is 
a small fraction of the mass of the star. For a light crust, 
17 <C 1 and g ^ GM/r^; for the incompressible case, 
U = 3 and g = GMr/R?. If g, is almost constant [or 
is replaced by an appropriately averaged value as below 


Eq. (68) of Ref. 0], Eq. 10 simplifies for an incompress¬ 
ible completely solid star to 

<522 ,max — 7^R®tTi„ax/(GM) (6) 

where 7 « 13. Evaluating Eq. 0 for a conventional neu¬ 
tron star with a thin crust and liquid core, 7 becomes 
about 120Ai?/i?, where Ai? « i?/10 is the thickness of 
the crust, and thus 7 is numerically almost identical. The 
appropriately averaged shear modulus from Ref. 0 is 
^ « 4 X 10^® erg/cm^, a factor of a few below its maxi¬ 
mum value at the bottom of the crust. 

Solid strange stars. —The idea that some “neutron 
stars” are in fact made of strange quarks was proposed in 
the 1970s 0. The idea that such stars are solid curre ntly 
is being pursued by Xu’s group, beginning with Ref. [l(l| . 
(This is distinct from a crystalline color superconduct¬ 
ing quark phase mi , which I do not consider here.) Xu 
notes that the burst oscillation frequencies observed in 
low-mass x-ray binaries correspond to the first few tor¬ 
sional modes of a solid strange star—if the matter has a 
typical shear modulus ^ « 4 x 10^^ erg/cm^, a thousand 
times the typical value in the crust of a conventional neu¬ 
tron star. Xu estimates that quarks clustered in groups 
of 18 or so could produce such a shear modulus. Since 
Ref. was published, the burst oscillation frequency 
has been observed to closely match the spin frequency of 
the neutron star in at least one system [ 13 . This renders 
the identification with torsional mode frequencies prob¬ 
lematic. However, the x-ray burst oscillation mechanism 
may be different for different binaries, and it is worth con¬ 
sidering the effect on the maximum elastic deformation 
if the shear modulus is very high for whatever reason. 

Using Xu’s shear modulus in Eq. © gives 



Bejger and Haensel find a different empirical formula 
for the moment of inertia for strange stars, 
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This combined with Eq. 0 yields a maximum ellipticity 
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for solid strange stars, where I have inserted the scalings 
of /i from Ref. M except for the / and the x depen¬ 
dence, which roughly cancel out. With self-gravity, the 
canonical number is Cmax = 6 x 10“^. 

Hybrid and meson condensate stars .—Glendenning 
m showed that the phase transition from baryonic mat¬ 
ter to quark matter occurs over a range of pressures 
rather than at a single value. (The argument holds for 
stars with charged meson condensates as well as for stars 
with quark-baryon cores 0 . The numbers are very simi¬ 
lar, so I discuss only hybrid stars.) Purely baryonic mat¬ 
ter at high densities is isospin asymmetric, which is ener¬ 
getically unfavorable. Moving toward isospin symmetry 
(creating more protons) would require negative charges 
to compensate, and leptons are not favored since they are 
nearly massless. When the quark phase becomes avail¬ 
able, baryonic matter can attain positive charge density 
by moving negative charge into areas of quark matter. 

The crystal structure of the mixed phase changes with 
density. Immediately above the threshold density for the 
beginning of the phase transition, the mixed phase con¬ 
sists of small quark droplets arranged in a bcc lattice in 
a baryonic background. As the density increases, the 
droplets grow and merge to become rods, then slabs. 
Eventually the baryonic matter becomes the minority 
slabs, then rods, then finally droplets before disappear¬ 
ing entirely. The locations of these layers are highly 
parameter-dependent; as an upper limit the mixed-phase 
crystal can occupy the innermost 8 km of the star . 

The shear modulus of a bcc lattice of point charges can 
be written in the parameters of Ref. Q as Hi 

= 0.075 q^D^/S'^. (10) 


I have assumed spherical droplets of (esu) charge density 
q, diameter D, and spacing S (Wigner-Seitz cell diam¬ 
eter, or '/3/2 times the lattice constant). Corrections 
due to the nonsphericity of the droplets reduce this by 
an amount that is small in most of the layer. Typical 
numbers from Chap. 9 of Ref. 0 give 

/r = 4 X 10^^ erg/cm^ 



This is of order 10^ times the typical value in the inner 
crust, mainly due to the charge of the droplets (about 10^ 
rather than Z < 55 in the crust) although the density of 
droplets is slightly greater too [ij . 

The dominant correction to Eq. d) is due to charge 
screening. This effect is difficult to evaluate precisely, 
but a rough estimate can be made as follows. Heiselberg, 
Pethick, and Staubo 0 estimate the screening length 


A = 




- 1/2 


( 12 ) 


in the mixed phase as 10 fm in the baryonic matter and 
5 fm in the quark matter. (Here Qi is the charge of 
species i, ni is its number density, and u, is its chemi¬ 
cal potential.) Detailed calculations [l^ of the partial 
derivatives in Eq. d for baryonic matter without a 
quark phase suggest that A ~ 5 fm is a lower limit. Since 
these lengths are comparable to the droplet size and sep¬ 
aration, screening will appreciably reduce electrostatic 
effects but not make them negligible. (The leptons can 
be neglected since their charge density is tiny y.) 

First note that screening does not appreciably change 
the droplet size. The quark volume fraction y « (D/S)^ 
is set by, e.g., the pressure, and D is found by minimizing 
the sum of surface and Coulomb energy densities at fixed 
X . The mean charge density q is then fixed even under 
rearrangement of charges due to screening. The Coufomb 
energy density can be written g{D/X)C{x)D"^, where g is 
a geometric factor. Coing from uniform density (D/A <C 
1) to a sheff of charge (D/A ^ oo) only reduces g by 
1/6. Using a rough approximation g oc 5 -I- exp(—D/A), 
the screened D is related to the unscreened Dq by 

-|-D^/(10A)exp(-D/A). (13) 

Since Dq and A are comparable, this is a few percent 
correction at most, and since (D — Dq) <C A, even the 
exponential factor is corrected only by a few percent. 

Civen that the droplet size does not change apprecia¬ 
bly, and given that the charges inside the droplet re¬ 
arrange themselves respecting spherical symmetry, the 
problem outside the droplet reduces to the classic screen¬ 
ing problem. The potential is multiplied by the Yukawa 
factor exp(—r/A). Since the shear modulus is roughly 
a second derivative of the potential energy, it is multi¬ 
plied by roughly this factor. For typical D and S values, 
screening then reduces the shear modulus by « 20 
(A « 10 fm) or e® ~ 400 for A « 5 fm. The up¬ 
per limit (weak screening) on the shear modulus is then 
/j, « 2 X 10®^ erg/cm®. 

The effective shear modulus for the rod and slab config¬ 
urations can be estimated from the droplet result. Matter 
made of rods cannot resist a shear stress along the axis 
of the rods, but will have a perpendicular response simi¬ 
lar to that of the droplets. This anisotropic case requires 
an elastic modulus tensor rather than a shear modulus 
scalar. However, if the glitch history of the neutron star 
has led to granulation [3j , the formation of small domains 
with different principal directions, then a macroscopic 
rms response of the matter averaged over many domains 
is isotropic with an effective shear modulus reduced by 
or y^ 1/3, which can be neglected here. 

Now evaluate Eq. ©• The density of the core of the 
star varies only by a factor of a few, so use the incom¬ 
pressible limit. Most of the integral comes from the 
droplet and rod layers, where the weak-screening shear 
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modulus is roughly constant at 2 x 10^^ erg/cm^. Then 
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where is the radius of the hybrid core. Bejger and 
Haensel |8| find that hybrid stars obey the same moment 
of inertia relation © as normal neutron stars, so 



for a fiducial value of 3 x 10“®, or up to 9 x 10“® with 
the self-gravity of the deformation. 

Implications .—What are the immediate consequences 
for LIGO? The S2 paper 0] quotes direct gravitational- 
wave observational upper limits on e for 28 pulsars. How¬ 
ever, 19 of these pulsars already have lower indirect upper 
limits on e (typically 10“® or less) due to the measured 
spin-downs. The remaining 9 pulsars are in globular clus¬ 
ters where the spin-down is obscured by acceleration, and 
thus have no competing upper limit. The S2 upper limits 
on e for these are 4-24x10“^, all within the maximum 1 
estimate for a solid strange star. With LIGO’s upcoming 
data run at full initial sensitivity, the same 9 pulsars will 
be observable at e of 1-8x10“®, within the maximum 
for hybrid stars; and the Grab pulsar will be observable 
at e = 1.2 X 10“"^, 6 times less than its spin-down limit 
and within the solid strange star range [ij. An all-sky 
search for unknown neutron stars could detect hybrid 
stars within a kpc and solid strange stars at the galactic 
core with tens of teraflops computing power 113 

If a pulsar is observed in gravitational waves with 
e ^ 10“^, it cannot be a conventional neutron star. 
An upper limit (nondetection) at higher e does not rule 
out any exotic model—a given star may happen to be 
nowhere near its breaking strain. However, with enough 
strict upper limits, population statistics and deformation 
mechanisms can be constrained. 

Bildsten 0 proposed that the spin frequencies of stars 
in low-mass x-ray binaries are set by equilibrium be¬ 
tween accretion torque and gravitational radiation from 
thermally induced deformations of the crust. In order 
to match the observed spin frequencies, this requires 
quadrupoles « 10®® g cm^. For the exotic stars consid¬ 
ered here, such quadrupoles under anisotropic accretion 
are possible if the breaking strain is smaller than 10“® 
or if the detailed accretion physics (temperature depen¬ 
dence, spreading of material, etc.) prevents achieving 
breaking strain. 

Gravitational waves from freely precessing neutron 
stars have been considered poor prospects even for ad¬ 
vanced LIGO. But if internal damping is weak, a popu¬ 
lation of stars precessing after birth with e = 10“^ would 


be detectable with broadband advanced LIGO 0| . 

A starquake that causes a glitch will also cause a burst 
of gravitational waves as the modes of the star are excited 
and ring down. This amplitude is determined by the 
energy in the glitch, which is determined by the observed 
frequency jump and thus is not affected by exotic matter 
models. But the maximum elastic energy in the star 
scales as the shear modulus, and thus is up to 10® times 
larger for quark stars than for conventional neutron stars. 
Vela glitches are still too large and frequent (by several 
orders of magnitude) to be explained as quakes, but the 
mean time predicted between quakes is reduced for the 
Grab pulsar to a few years |3| —comparable to what is 
observed. 

After this Letter was submitted, SGR 1806-20 under¬ 
went a giant superflare 0113 with estimated energy 
more than 10"^® erg. Theoretical models equate this en¬ 
ergy with the maximum elastic energy of the star, which 
is problematic for a normal crust but feasible with the 
exotic models considered here. 
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